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SUMMARY 

A computational method/procedure Is described which can be used to simu- 
late Individual and mixed mode Interlaminar fracture progression In fiber com- 
posite laminates. Different combinations of Modes I, II, and III fracture are 
simulated by varying the crack location through the specimen thickness and by 
selecting appropriate unsymmetrlc laminate configurations. The contribution of 
each fracture mode to strain energy release rate Is determined by the local 
crack closure methods while the mixed mode Is determined by global variables. 
The strain energy release rates are plotted versus extending crack length, 
where slow crack growth, stable crack growth, and rapid crack growth regions 
are easily Identified. Graphical results are presented to Illustrate the 
effectiveness and versatility of the computational simulation for (1) evalu- 
ating mixed-mode Interlaminar fracture, (2) for Identifying respective dominant 
parameters, and (3) for selecting possible simple test methods. 


INTRODUCTION 

Interlaminar delamination of angleplled laminates Is a fracture mode which 
needs to be carefully examined and properly accounted for In the design of com- 
posite structures. Regions prone to delamlnatlons Include free edges, loca- 
tions of stress concentration, joints, Inadvertently damaged areas, and defects 
arising during the fabrication. One way to properly account for Interlaminar 
delamination In design is to determine Interlaminar fracture toughness para- 
meters and then to compare these parameters to their respective critical values 
and to stress states which are likely to Induce Interlaminar fracture. 

Interlaminar fracture In angleplled laminates Is generally Induced by 
Individual and/or mixed mode type (Modes I-openlng, II-shear1ng, and Ill- 
tearing) fracture. In order to properly assess composite fracture resistance, 
the fracture toughness parameters for each mode and for mixed modes must be 
known or determined. These parameters can be determined either experimentally 
or computationally. The unsymmetrlc double cantilever, the mixed mode flexure 
and end-notch flexure specimens can be used to experimentally measure mixed 
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mode fracture toughness. Variations of these test methods can be used to indi- 
vidually measure Mode I and Mode II. Mode III, on the other hand, is usually 
measured as a combination with either Mode I or Mode II since Mode III is gen- 
erally a very difficult test to perform. 

A computational procedure (ref.l) was developed at NASA Lewis Research 
Center for predicting interlaminar fracture in unidirectional composites. This 
computational procedure consists of three-dimensional finite element analysis 
in conjunction with composite micromechanics. It is used to determine fracture 
toughness parameters by computationally simulating respective tests as follows: 

Double Cantilever Mode I 

End-notch-flexure Mode II 

Mixed-mode Flexure Mixed Mode I and II 

Recently this computational procedure has been modified to determine mixed 
Modes I, II, and III in composite laminates. The objective of this report is 
to describe this modified procedure and its application to composite inter- 
laminar mixed mode fracture. 

The modified procedure consists of three-dimensional finite element anal- 
ysis in conjunction with Integrated composite mechanics (micromechanics, 
macromechanics, combined-stress failure criteria, and laminate theory). The 
procedure is used to computationally simulate the mixed-mode fracture of flex- 
ural specimens made from unsymmetrlc/unbalanced laminate configurations. Dif- 
ferent combinations of Modes I, II, and III are simulated by varying the crack 
location through the specimen thickness and by selecting appropriate unsym- 
metric laminate configurations. The contribution of each fracture mode Is 
determined by local crack closure methods while the mixed mode is determined 
from the global method. The fracture modes are determined in terms of strain 
energy release rates. The strain energy release rates are plotted versus 
extending crack length, where slow crack growth, stable crack growth, and rapid 
crack growth regions are easily identified. Graphical results are also pre- 
sented which illustrate the effects of parameters such as: (1) ply orienta- 

tion, (2) laminate configuration, (3) interlaminar crack location, and (4) 
laminate material-coupling coefficients on strain energy release rates. Fund- 
amental considerations and possible generalizations are also described. 


COMPUTATIONAL SIMULATION: FUNDAMENTAL CONSIDERATIONS 

The computational simulation method for evaluating interlaminar individual 
and mixed mode fracture toughness consists of three-dimensional finite element 
analysis, including finite element local mesh refinement, and integrated com- 
posite mechanics. Several fundamental considerations are associated with this 
computational simulation method. These fundamental considerations include: 

(1) laminate configurations, (2) component geometry and loading, (3) finite 
element model, (4) composite system, and (5) computational procedure. Each of 
these are described below. 


Laminate Configurations 

The laminate configurations used in these studies were unbalanced, 
unsymmetric These laminates were selected in order to evaluate 
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the effects of the different laminate material-coupling coefficients. In the 
laminate force deformation relationships, on the Individual and mixed mode 
fracture strain energy release rates (SERR). 

The force deformation relationships are given by the following matrix 
equation (ref. 2 ): 
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The notation In equation (1) Is as follows: 

{N} denotes the section (know or calculated from eq. (1)) membrane or axial 
forces (N xx , Nyy and N xy ) along the respective structural axes; {M} denotes 
the corresponding bending moments; [A], [C], and [D] denote membrane, coupled 
bending-membrane and bending stiffness [3 x 3] matrices, respectively; {e 0 } 
denotes the mid plane strains U xxo » c yy 0 » c xyo ); *** denotes the laminate 

structural axes curvatures (k xx , K y y, and K xy ); AM and AT denote changes 
In moisture and temperature, respectively; {Nm> and (Mm) denote axial 
forces and bending moments due to moisture gradient through-the-lamlnate- 
thlckness; and {Ny} and {My} denote forces and moments due to correspond- 
ing temperature gradient. The material-coupling coefficients which couple the 
different fracture modes are: (1) for N xx /N xy (Mode II/Mode III) cou- 
pling, (2) A 23 for Nyy/N xy (Mode II/Mode III) coupling, (3) C 13 for 

Nxx/Mxy (Mode II/Mode I and Mode III) coupling, ( 4 ) C 23 for Nyy/M X y (Mode 
II/Mode I and Mode III) coupling, (5) D *| 3 for M xx /M X y (Mode I and Mode II/ 
Mode II and Mode III) coupling, and ( 6 ) D 23 for M yy /M xy (Mode I and Mode II/ 
Mode II and Mode III) coupling. 

The specific laminates studied are summarized In table 1. Seven cases 
( t _e 36^ f ®12^ » 0 = 0°» 15°, 3 0 o , 45°* 60°, 75°, and 90°) were Investigated 
for the laminate configuration effects on the SERR. At the © = 45° position, 
a total of three cases ( [ -45 m /45y 2 ] * m * 36, 60 and 84) were Investigated 
for the Interlaminar crack location effects on SERR. It Is Important to note 
that these specific laminate configurations and the Interlaminar crack loca- 
tions were selected only for computational simulation convenience. They repre- 
sent just one application of the present procudure. 


Component Geometry and Loading 

The geometric configuration and the loading of the component selected for 
the studies are shown In figure 1. This component, dimensions and loading were 
selected because: ( 1 ) It Is the same specimen used for measuring Interlaminar 

and mixed mode (I and II) fracture toughness as described In reference 1; (2) 

It Is one of the simplest component-loading condition combinations that can be 
used to computationally evaluate the effects of laminate configuration, 
material-coupling and Interlaminar crack (delamlnatlon) location on Individual 
(Mode II) and mixed mode fracture (Modes I and II and Modes I, II, and III) 
SERR; and (3) It Is readily amenable to conduct experimental studies for com- 
posite system screening and the determination of environmental effects. 
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Finite Element Model 


The computational simulation was performed using MSC/NASTRAN with local 
mesh refinement. A computer plot of the finite element model Is shown In 
figure 2. The entire component with the supports and typical overhang Is 
modeled. The finite element model consisted of 1856 solid elements, 2450 
nodes, and 7350 degrees-of-freedom (DOF). This finite element model Is simi- 
lar to that used In reference 1 but with variable thickness elements. A fron- 
tal view of the finite element model showing the local details and the location 
of the mesh refinement "super element" Is shown In figure 3. A schematic of the 
local superelement Is shown In figure 4 where Its characteristics are also sum- 
marized. It Is worth noting that the local model Includes solid elements for 
the Interlaminar matrix layer. The Inclusion of the Interlaminar matrix layer 
and the anaylsls of the entire specimen are two unique features which distin- 
guish the computational simulation method described herein from what Is conven- 
tionally done and reported In the literature. 


Composite System 

The laminates are assumed to be made from a composite material consisting 
of AS graphite-fibers In an Intermediate-modulus high-strength epoxy (AS/E) 
with a 0.6 fiber volume ratio (FVR). The composite material properties 
required for the three-dimensional finite element analysis were generated 
using the (Integrated composite analyzer) computer code (ICAN) rising typical 
(ref. 2) constituent material properties summarized In table 2. These proper- 
ties are encoded as subsets of the ICAN resident data bank. The three- 
dimensional MSC/NASTRAN material card properties generated by ICAN are summa- 
rized In table 3 using the NASTRAN notation. Note that the last entry In 
table 3 Is the interlaminar (Interply) layer thickness which Is predicted by 
ICAN and which Is the finite element thickness for this layer. Using Inte- 
grated composite mechanics to predict the requisite three-dimensional finite 
element material properties and the Interply layer thickness are two additional 
unique features of the present computational simulation method. 


Computational Procedure 

The computational procedure used Is that developed previously (ref.l) 
modified to accommodate unsymmetrlc laminates for Mode III fracture. It con- 
sists of a global method for determining mixed mode fracture SERR and the local 
crack closure method for determining the contributions of the Individual modes 
fracture SERR. Each method Is summarized below for completeness. 


Global Method 


The specific computational steps for this method are as follows (refer to 
fig. 1): 

(1) Model the component with crack length "a" using three-dimensional 
finite elements. 

(2) Apply a load (P) at component midspan. 
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(3) Calculate the midspan displacement w(a) (w as a function of "a") 
using three-dimensional finite element analysis (FEA). 

(4) Induce crack extension Aa keeping load (P) constant. 

(5) Calculate the midspan deflection [w(a ♦ Aa)]. 

(6) Determine the Strain Energy Release Rate (SERR), Gy, from 

c - p [w( a ♦ . Aa) - w(a ).] 
b T ‘ y 2b Aa 

where b Is the specimen width. 

(7) Repeat steps (4) to (6). 

(8) Plot results for Gy versus a or Aa. 

(9) Identify fracture toughness characteristics as described In 
reference 1 . 

(10) Examine complete stress state near crack tip. 

(11) Compare with corresponding uniaxial composite strengths. 

(12) Look for possible correlation of fracture toughness parameters with 
composite uniaxial strengths. 

The global method yields the global "fracture toughness" SERR without any 
regard to participating and/or dominating local fracture modes. 

Local Closure Method 

The specific steps for this method are as follows: 

(1) Perform steps (1) and (2) as In the Global Method. 

(2) Calculate displacements (u(a), v(a), w(a) at the crack tip nodes 
using FEA. 

(3) Induce crack extension Aa keeping P constant. 

(4) Calculate (u(a + Aa), v(a + Aa), w(a + Aa)) at the same nodes as In 
step (2). 

(5) Apply enforced displacements (single point constraints) using the 
step (2) displacements (u(a), v(a), w(a)) at the crack tip nodes. 

(6) Repeat FEA with these single point constraints. 

(7) Calculate the corresponding forces at these constraints (F x , Fy, 

F z ). These are called the single point constraint forces In FEA. 
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(8) Determine the local SERR's from 


G . - F z 


fwU * Aa) - w(a)l 
2b Aa 


'II 


r ru(a + &a) - u(a)l 
x 2b Aa 


(3) 

(4) 

(5) 


G _ F [ v(a + Aa) - v(a) ] 
b III ‘ r y 2b Aa 

(9) Repeat steps (3) to (8). 

(10) Follow steps (8) to (12) In the Global Method. 

The local crack closure method yields the contribution of each local frac- 
ture mode to the composite mixed mode fracture toughness (SERR). 


RESULTS AND DISCUSSION 

The primary result of this study Is the demonstration of the effectiveness 
and versatility of the computational simulation method. The results presented 
herein evaluate composite structure and material factors that Influence 
Individual and mixed mode Interlaminar fracture. These factors Include: 

(1) Interlaminar crack opening (a), (2) ply orientation (e), (3) laminate 
configuration ([-©/+©]), (4) Interlaminar (Interply) crack location ([©„,/+©„]), 
(5) material-coupling coefficients, and (6) three-dimensional stress state 
ahead of the crack tip. Each of these factors are described below In terms of 
their effects on Individual and mixed mode fracture SERR relative to component 
and loading In figure 1. 


Effects of Crack Opening 

The effects of crack opening on the Individual mode fracture SERR (Gj, G j j 
and Gin) and mixed mode SERR (Gt) are plotted In figure 5 for the seven dif- 
ferent laminate configurations, [-©36/ + ©12]. ( e = 0°. 15°, 30°, 45°, 60°, 75° 
and 90°), figures 5(b) to (f), respectively, assuming a 480 lb load. The 
Important observations In figures 5(b) to (f) are: (1) the shearing fracture 
Mode II ( Gj j ) dominates the stable crack propagation (growth); (2) the open- 
ing fracture Mode I ( Gj ) dominates the rapid crack propagation; (3) the tear- 
ing fracture Mode III (Gm) has generally similar form as the shearing mode 
(Gji) but of considerably lower magnitude (about 10 percent), and (4) the mixed 
mode fracture (Gj) Is not the algebraic sum of the other three and, therefore, 
needs to be determined using a global method. 

The Important conclusions are: (1) a global method should be used to 

determine the mixed mode fracture, (2) local methods effectively determine the 
contribution of the Individual fracture modes, and (3) local methods can be 
used to determine the crack opening range In which the different fracture modes 
dominate. One significant Implication from the above observatlons/concluslons 
Is that local averaging methods, such as: crack opening, J Integral, stress 
Intensity, Intense energy parameters, and Inherent material flaw parameters 
used to estimate global fracture parameters, will generally underestimate the 
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mixed mode fracture as defined by SERR Gy. Use of these local averaging 
methods will result In an optimistic estimate of the composite structure Inter- 
laminar fracture toughness. The loss In energy In the removed Interply layer 
elements contributes about three percent to this optimistic estimate. Experi- 
mental data obtained by using local measuring techniques should be Interpreted 
with the above Implication In mind. 


Effects of Ply Orientation 

The effects of ply orientation on the maximum Individual and mixed mode 
fracture SERR (assuming a 480 lb load) are plotted In figure 6 for the seven 
cases of the [-© 36 A© 12 ] (© = 0% 15°, 30°, 45°, 60°, 75° and 90°) AS/E 
laminate. The relative dominance of the opening fracture mode SERR (Gi) and 
the negligible contribution of the tearing fracture mode SERR (Gm) on the 
maximum mixed fracture mode SERR (Gy) are clearly observed In this figure. 
Another observation Is that the "maximum" magnitude of mixed mode fracture (Gy) 
levels off at ply angle orientations greater than 60° (0 > 60°). 

The above observations lead to the following conclusions: (1) the rapid or 
unstable Interlaminar crack growth Is dominated by the opening fracture mode; 
(2) the tearing fracture mode SERR (Gjij) Is negligible compared to Gy for ply 
orientation angles greater than 60°; (3) the Individual fracture modes (Gj and 
Gn) and the mixed mode fracture are practically Independent of ply orientation 
angle greater than 60°; and (4) ply orientation angles less than 60° have sig- 
nificant Influence on the SERR of Individual fracture mode and to mixed mode 
fracture. 


Effects of Laminate Configuration 

The effects of the laminate configuration on the SERR due to a 480 lb load 
are shown In figure 7. The effects on the Individual fracture mode SERR's are 
plotted In figures 7(a) to (c). The SERR for the mixed fracture mode (Gy) Is 
plotted In figure 7(d). The Important observations In this figure are: (1) 

the shearing mode fracture SERR ( Gy j ) appears to reach a maximum and then 
decrease with Increasing crack opening for practically all the laminate config- 
urations, (2) the tearing mode fracture SERR (Gm), on the other hand, con- 
tinues to Increase with crack length for some laminate configurations. 

The conclusion from the above discussion Is that laminate configurations 
can be selected for "stable" shearing and tearing fracture mode crack growth 
for given composite components and loadings. It Is Important to keep In mind 
that advantages of this can be taken only In the absence of opening mode 
fracture. 


Effects of Interlaminar Crack Location 

Recall that the effects of the Interlaminar crack location are studied 
using the laminate [-©^/+©t 2 ] . For these studies © Is chosen to be 45° 

since the laminate configuration with this ply angle had the greatest magni- 
tude of tearing mode fracture SERR ( Gy j j ) (see fig. 7(c)). Results obtained 

are shown In figure 8 for three different crack locations (m = 36, 60, and 84) 
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as defined In table 1. The loads corresponding to these crack locations were: 
480, 882, and 1358 lb, respectively. In order to subject the laminates to 
approximately constant external work. 

The effects of crack location on the Individual mode fracture SERR (Gj, 

Gjj and Gni) are plotted In figures 8 (a) to (c). That for the mixed mode 

fracture SERR (Gy) Is plotted In figure 8 (d). Collectively these results show 
that the effects of the Interlaminar crack location on the SERR decrease as the 
thickness above the crack plane Increases. These effects are more pronounced 

for the opening and tearing modes fracture SERR than they are for the shearing 

mode fracture. 

The above discussion leads to the following Important conclusion: lami- 

nates subjected to bending can sustain relatively large Interlaminar cracks 
when these cracks are located near the tensile surface of laminate. The sig- 
nificant Implication of this conclusion Is that relatively large Interlaminar 
crack sizes are detectable by available NDE Inspection methods. And, there- 
fore, this computational simulation methods may be used to establish guidelines 
for setting fracture control requirements, for example, size of allowable 
delamlnatlon versus depth from surface. 


Effects of Material-Coupling Coefficients 

The tearing mode fracture Is present In unbalanced, unsymmetrlc laminates 
when subjected to bending loads as was already mentioned. The magnitude of the 
tearing mode fracture SERR depends strongly on the magnitude of material- 
coupling coefficients such as A 23 , and C 23 . However, these coupling coeffi- 
cients generally are present In combinations with their complements and with 
other material-coupling coefficients A-j 3 , C-j 3 , and D 23 . 

The effects of ply orientation In a [-036/+©i2l laminate on the material- 
coupling coefficients are shown In figure 9. The effects on axial normal-shear 
coupling A -|3 and A 23 are plotted In figure 9(a). The effects on membrane- 
bending coupling Ci 3 and C 23 are plotted In figure 9(b). The effects on 
bending-twisting coupling Di 3 and O 23 are plotted In figure 9(c). As can 
be seen In figure 9 the 13 and 23 coupling coefficients are complementary, and 
vanish at 0 = 30° or 90°. 

The effects (variations) of the coupling coefficients on the opening mode 
fracture SERR (Gj) due to a 480 lb load are shown In figure 10. The effects 
of A 13 , Ci 3 and D 13 are plotted In figure 10(a), (b), and (c), respec- 
tively. As can be seen In these plots Gi Is a double-value function of the 
material-coupling coefficients. This double-value function results from the 
complimentary parts of the material-coupling coefficients mentioned previously. 

Comparable plots for the effects on the shearing mode fracture SERR ( Gj j ) 

are shown In figure 11. These effects are similar to those for the opening 
mode. The effects of the material-coupling coefficients on the tearing mode 
fracture SERR (Gm) are shown In figure 12. These effects are also double- 
value functions and reach their maximum magnitude at 0 = 45°. 

The Important conclusion from the results In figures 11 and 12 Is that 
laminate configurations exist which exhibit considerable tearing mode fracture 
SERR (Gih) flexural components which are subjected to bending loads. It 
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Is possible, therefore, to determine the magnitude of this mode by simple 
experimental techniques. This magnitude will not be measured Individually. It 
will be In combination with the other two modes. 


Three-Dimensional Stress State 

The Individual fracture modes are caused by Individual, or combination of 
stresses at the crack tip. The dominant stress for each respective mode can be 
determined by plotting the three-dimensional stress state versus distance from 
the crack tip. The three-dimensional stress state Is computed as a part of the 
three-dimensional finite element analysis. This procedure Is Illustrated In 
figure 13 where the three-dimensional stress state In the Interply layer ahead 
of the crack tip Is shown for the [ -30 36 /-«-30 i 2 ] laminate. 

The Important observations In figure 13 are: (1) o zx has the highest 

magnitude (about -25 ks 1 ) ; (2) o zz has the second highest magnitude (about 
14 ks 1 ) ; the magnitudes of the other stresses In decreasing order are as 
follows: o xx approximately 8 ksl, o yy approximately 5 ksl, o yz approximately 

4 ksl and o xy approximately 0. The dominant stresses, therefore, which cause 
Interlaminar fracture growth and their corresponding fracture modes are as 
follows (dominant-stress/fracture-mode) : 

1. "Longitudinal" Interlaminar shear stress/shearing fracture mode 
^zx /G Il) 

2. Interlaminar normal stress/opening fracture mode 
(<j zz /Gi) 

3. "Transverse" Interlaminar shear stress/tearing fracture mode 

(oy Z /Gin) 

The In-plane stresses also contribute to the Individual fracture modes. 
Following the above notation their contributions are: 


1. 

^xx^II 

and 

G I 

2. 

<Jyy/GlH 



3. 

<J X y = 0 

utes to 

and. 

Gill 

therefore, does not contribute; otherwise It contrlb 


The plots In figure 13 can also be used to estimate stress magnitudes at 
which crack propagation will occur. This Is accomplished by plotting the cor- 
responding ply strengths. For example the longltudlng Interlaminar shear 
strength for this composite Is about 13 ksl. The transverse Interlaminar shear 
strength Is about 9 ksl. The normal Interlaminar stress Is about 6 ksl. Com- 
paring the dominant stress magnitudes It Is found that: 

o zx approximately 25 ksl > 13 ksl; a zz approximately 14 ksl >7 ksl; and 
oy Z approximately 5 ksl < 9 ksl. Therefore, It may be concluded that crack 
propagations will occur due to the shearing and opening fracture modes. 
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SUMMARY 


The significant results and conclusions of an Investigation to computa- 
tionally simulate composite Interlaminar Individual and mixed-mode fracture 
toughness as determined by strain energy release rates (SERR) are summarized 
below. 

1. The Individual and mixed mode fracture SERR can be readily determined 
using a computational simulation procedure that consists of three-dimensional 
finite element analysis and Integrated composite mechanics. 

2. Individual and mixed mode fracture SERR magnitude of [-0^/+©,,] are 
strongly Influenced by crack length, ply angle, and Interlaminar crack loca- 
tion. However, the maximum magnitude of the mixed mode fracture SERR (Gy) 

Is practically Independent for ply orientations greater than 60°. 

3. The tearing fracture mode SERR (Gju) has the smallest magnitude compared 
to opening (Gi) and shearing (Gn) fracture modes for this case. The tear- 
ing fracture mode Is generally present In combinations with other fracture 
modes . 

4. The Individual and mixed mode fracture SERR decrease as the location of the 
Interlaminar crack approaches the tensile surface of a laminate which Is sub- 
jected to bending. 

5. The magnitudes of the material coupling coefficients strongly Influence the 
Individual and mixed mode fracture SERR. The presence of these coefficient 
magnitudes Induce tearing mode fracture In combination with other fracture 
modes when the component or specimen Is subjected to In-plane or bending loads. 

6. The dominant stresses associated with Individual fracture modes are readily 
Identified from the three-dimensional finite element analysis results. 

7. Laminate configurations can be selected which exhibit substantial tearing 
fracture mode In a three-point-bend specimen, and therefore, the magnitude of 
this mode can be determined by simple experiments. 

8. Stress magnitudes ahead of the crack tip can be compared with corresponding 
local laminate strengths In order to determine the dominant stress which drives 
the crack. 

9. Local averaging methods/techniques may provide optimistic estimates of 
global fracture toughness critical parameters. 

10. Collectively the results demonstrate that the procedure described herein 
can be used to computationally simulate/evaluate mixed mode fracture tough- 
ness parameters In composite components subjected to complex loadings. 
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TABLE I. - LAMINATE CONFIGURATIONS 


No. 

Configuration 

Total 

thickness , 
In. 

Crack location 
from bottom 
surface, 

In . 

1 . 

‘W 

0.24 

0. 

06 

2. 

t- 1 5 36 A15 12) 





3. 

[-30 36 /. 3 0i 2 ] 





4. 

‘-♦ V 45 !* 1 





5. 

t-60 36 /.60 12 ] 





6. 






7. 

'- 9 V t 9 V ] 





8. 

[-45 60 A«5 12 ] 

0 . 

36 

\ 


9 . 

C- 4 V t45 12 ] 

0. 

48 

0. 

06 


TABLE II. - TYPICAL CONSTITUENT MATERIAL PROPERTIES USED IN THE 
SIMULATION 


Property 

Units 

Fiber 

AS-graphlte 

Matrix 

Intermediate-modulus 

high-strength 



Symbol 

Value 

Symbol 

Value 

Elastic moduli 

xlO 6 psl 

Efll 

31 

Em 

0.5 



Ef 22 

2 





Of 12 1 

2 





Gf23 

1 

Gm 

.185 

Poisson's ratio 


Vf12 

0.2 

v m 

.35 



Vf 23 

0.25 



Strengths 

ks 1 

SfT 

400 

SmT 

15 



$fc 

400 

$mc 

35 

Fiber diameter 

Inch 

<*f 

0.0003 

S ms 

13 





TABLE III. - MATERIAL PROPERTIES FOR THREE-DIMENSIONAL FINITE ELEMENT ANALYSIS 
[ICAN predicted, reference 2.] 


Property 

Units 

Symbol 

(MSC/NASTRAN) 

Laminate configuration 

0 

±15 

±30 

±45 

±60 

±75 

±90 

Elastic 

rapsl 

6 n 

19.080 

16.840 

11.500 

6.035 

2.692 

1.589 

1.469 

constants In 


6 21 = 6 12 

0.546 

1 .607 

3.728 

4.788 

3.72P 

T .607 

0.546 

the stress/ 


G - G 

31 13 

0.546 

0.552 

0.567 

0.589 

0.610 

0.626 

0.631 

strain 


G 41 = °14 

0.000 

±4.039 

±5.650 

±4.404 

±1.968 

±0.365 

0.000 

relationship 


G 51 ’ G 1S 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

for AS/E 


G 61 - G U 

0.000 

0.000 

0.000 

0.000 
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Figure 2.- Computer plot of component 3-D finite element model (1856 solid elements 
32-Pentahedrons, 1824 Hexahedrons, 2450 nodes, 7350 degrees-of-freedom (DOF)). 
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Figure 3.- Frontal view of finite element model show- 
ing LOCATION AND GRID OF SUPERELEMENT. 



Figure 4.- Enlarged frontal view of superelement 
(360 SOLID ELEMENTS: 32-Pentahedrons, 328- 
Hexahedronj 450 NODES, 135Q DOF). 
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Figure 5.- Effects of crack opening on 

STRAIN ENERGY RELEASE RATES (INDIVIDUAL 
AND MIXED MODES) FOR DIFFERENT LAMINATE 
CONFIGURATIONS : AS“GRAPHITE 

FIBER/EPOXY (AS"EPOXY), WITH 0,6 FIBER 
VOLUME RATIO (FVR>) (A) -0 = 0°; 

(B) -0 = 15°; (c) -0 = 30°; (D) -0 = ^5°; 
(E) -0 - 60°; (F) -0 = 75°; (G) -0 = 90°. 
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Figure 5. - Concluded. 
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IGURE 6.- EFFECT OF PLY ORIENTATION ON 
MAXIMUM INDIVIDUAL AND MIXED MODE STRAIN 


ENERGY RELEASE RATES ([-0vf-/+0 
EPOXY WITH 0.6 hVR) . 


12 J; AS/ 
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Figure 7.- Effect of Laminate configuration 
on strain energy release rates <[-03g/ 

+0^1 2 1 • AS/E with 0.6 FVR) . a-Opening 
MODE (Gj) B-SHEARING MODE (Gjj), C-TEARING 
MODE (G hi ), D-MlXED MODE (Gj). 




D13 AND D23, LB-IN. C12 AND C23. LB A13 AND A23. LB/IN. 





Figure 9 .- Effect of ply orientation on 
Laminate material -coupling coefficient 

([_0 36 /+0 12^‘ AS/t WITH 0,6 hVR) - 
a-Axial normal-shear coupling A 13 and A 2 3 . 

B-MEMBRANE-BENDING COUPLING C 13 AND C23; 
C-BENDING-TWISTING COUPLING AND D23. 
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Figure 11.- Effects of Laminate material- 
coupling COEFFICIENTS ON SHEARING MODE 
STRAIN ENERGY RELEASE RATE (Gjj), 

([ ’ 0 36 /+e l2 ] ' AS/E WITH °' 6 FVR)< 
A-COUPLING COEFFICIENT (LB/IN.); 

B-COUPLINE COEFFICIENT (LB); 

C-COUPLING COEFFICIENT (LB-IN.). 






h I CURE 13.- 3-D STRESS STATE IN THE INTERLAMINAR MATRIX LAYER AHEAD OF CRACK TIP 
<[-30 36 /+30 12 ]. AS/E with 0.6 FVR). 
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